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ABSTRACT 

Wc study the effect of entropy injection in the intracluster medium (ICM) in light of 
the recent observationally determined universal pressure profile of the ICM. Beginning 
with a power-law entropy profile that is expected in the absence of any feedback, we 
show that a simple universal prescription of entropy injection results in the final, 
observed universal pressure profile. This simple prescription has two components, one 
associated with an overall increase in entropy and another associated with injection in 
the central parts of the cluster. Importantly, both the components of entropy injection 
are needed to produce the final universal pressure profile. This is indicative of a need 
of both preheating the ICM as well in situ AGN/SNe heating. We demonstrate the 
usefulness of the method by extending the calculations to clusters at high redshift, and 
predict redshift evolution of cluster scaling relations that can be tested against data. 
We show that the self-similar evolution of the universal pressure profile is equivalent 
to a negative evolution of entropy injection with redshift, with a scaling Sinj oc (1 -|- 
z)~°-^Sinj{z = 0). We also show the current observational data is indicative of the 
entropy injection decreasing with redshift. 

Key words: galaxies: cluster: general ~ galaxies: intergalactic medium - X- 
rays: galaxies: clusters - cosmology: miscellaneous 



1 INTRODUCTION 

Galaxy clusters are being increasingly used as cosmologi- 
cal probes. The evolution of galaxy cluster abundances, for 
example, is sensitive to the values of cosmological parame- 
ters, and have been extensively studied in order to determine 
or constrain these parameters (Wang & Steinhardt 1998, 
Holder et al. 2001, Majumdar & Mohr 2003, 2004, Gladders 
et al. 2007, Vikhlinin et al. 2009, Rapetti et al. 2009, 2010, 
Rozo et al. 2010, Sartoris et al. 2010). These attempts have 
been, however, limited by the inadequacy of our knowledge 
of the state of the intracluster medium (ICM) of clusters. 
Cluster scaling relations show that the distribution of the 
ICM gas is different from simple models of cluster forma- 
tion in which baryonic gas only reacts to the gravitational 
potential of the dark matter (Benson et al. 2004, da Silva et 
al. 2004, Bonaldi et al. 2007, Bonamente et al. 2008). 

It is now believed that a number of non-gravitational 
processes may have affected the ICM gas, including radiative 
cooling and energetic feedback from active galactic nuclei 
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(AGN) and/or supernovae (for example, see Voit & Bryan 
2001, Nath & Roychowdhury 2001, Balogh et al. 2001, Bor- 
gani et al. 2004, Roychowdhury et al. 2004, Kravtsov et al. 
2005, Puchwein et al. 2008, Battaglia et al. 2010). These pro- 
cesses tend to raise the entropy of the ICM gas, essentially 
heating the gas and pushing it outward, and they make the 
gas less X-ray bright for a given temperature, and help in 
reconciling theoretical predictions with observed data. 

In recent years, the Sunyaev-Zel'dovich (SZ) tempera- 
ture distortion of the cosmic microwave background radia- 
tion (CMBR) has received a lot of attention, in light of the 
possibility of observing a large number of clusters through 
SZ effect, since the magnitude of the effect is independent 
of redshift (Birkinshaw 1999). The scaling relations for SZ 
effect is also altered by the non-gravitational processes men- 
tioned above, and in general the magnitude of the SZ effect 
is decreased by the non-gravitational processes mentioned 
above (White, Hernquist & Springel 2002, Roychowdhury, 
Ruszkowski & Nath 2005). 

Recently, a universal pressure profile for the ICM gas 
has been shown to fit the X-ray observations of clusters (Ar- 
naud et al. 2010). This pressure profile can be combined with 
a mass model for the dark matter halo to yield the corre- 
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spending gas density and temperature profiles, and therefore 
its entropy profile, fn tfiis paper, we compare tfiis entropy 
profile of tfie ICM (inferred from tlie observed universal pres- 
sure profile) to determine the additional entropy needed for 
the ICM gas, starting with a physically motivated profile 
that is expected in the absence of any feedback. The aim is 
to determine the additional entropy needed to explain the 
observations of local clusters in a phenomenological manner. 

We then speculate upon the possible evolution of en- 
tropy injection with redshift, and calculate the correspond- 
ing pressure profiles at high redshift. These profiles can be 
tested against the data for the evolution of cluster scaling 
relations with redshift (e.g, between X-ray luminosity and 
temperature), and the evolution of entropy injection with 
redshift can be inferred with the help of this method. We 
show that although present available data are insufficient to 
constrain the evolution of entropy injection, one will be able 
to iterate with better data in future and determine the en- 
tropy injection at high redshift from observed cluster scaling 
relation evolution. 

The rest of the paper is structured as follows: In §2, we 
show different ways to calculate the ICM pressure profiles, 
in §3 we model injected entropy, in §4 we look at its evolu- 
tion with redshift and we finally conclude in §5. Throughout 
the paper, we use the subscript U for physical parameters 
related to the universal pressure profile by Arnaud et al. 
(2010). Also, we use the subscripts i and / for initial and 
final profiles and values. 



2 ICM PRESSURE PROFILE 

Recent observations by Arnaud et al. (2010) have shown 
that a scaled universal profile for the pressure of the ICM 
can fit the X-ray observations. Defining x = r/7?5oo, where 
-Rsoo is the radius within which the mean over density is 500 
with respect to the critical density, this profile is given by 
(eqns 9, 11, 13 of Arnaud et al. 2010), 



Pu{x) = Pst 
where, 

P50o(keV cm" 



Po hlr 
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(1) 
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3 X 1014 hg l Me 
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8.403 hg'^^^ , 
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where C500 



1.177, a 



1.051, P = 5.4905, 7 = 0.3801, 



ap = 0.12, a'p{x) ^ 0.1 - {a^ + 0.1) j^^^, M500 is the 
mass enclosed within -R500, H — 70 ho,7 km s^^ Mpc^^, and 
E(z) = H{z)/Ho = [^0(14- z)^ + Ao]^''^ is the ratio of the 
Hubble constant at redshift z to its present value. 

We assume a Navarro-Frenk- White (NFW) profile for 
the dark matter distribution (Navarro et al. 1997), along 
with an estimation of the concentration parameter 'c' from 
observations by Comerford & Natarajan (2007), given by 
c„,, = iMi^(Mvi./1.3 X 10^3 h-^ A/q)-"-i^±°-i«. We note 
that this prescription refers to the total cluster mass and 
is an empirical one. However, we have found that our fi- 
nal results are largely insensitive to the assumptions of the 



concentration parameter (see also Chaudhuri & Majumdar 
2011). 

We also use the analytical fit by Bryan & Norman 
(1988) for the over density Ac{z) = 18-k^ + 82x - 39 , 
with X — Qrn{z) — 1, with which one can calculate the virial 
radius, as. 



Mvir 



{4n /3)pcrit{z) Ac{z) 



1/3 



(3) 



In order to determine the additional entropy we need 
to use a default profile for the ICM gas that is predicted 
in the absence of any energy input. Following the tradition 
in X-ray astronomy, we define the 'entropy', S as kBT/nl^'^ 
(which is usually expressed in the units of keV cm^). 

There are three possibilities in this regard: 

• (a) One can either use an entropy profile that one ex- 
pects from the consideration of ICM gas getting heated at 
the accretion shock, and in the absence of radiative cooling. 
Voit et al. (2005) showed that simulations with only gravi- 
tational processes yield a power law entropy profile of type 
S,{r) = 1.32 X S200 (r/P2oo)' \ where 



S200 = 



2/3 5/3 



27r G'Ma, 



2/3 



(4) 



.15 hH{z) 

Here, we use an equivalent expression that uses P500, taken 
from Pratt et al. (2010), using the observed mean value of 
C500 = 3.2 for relaxed clusters (Pointecouteau et al. 2005), 



S{r)/S5oo = 1.42(r/i?5oo)''' 



where 



5*500 = lOekeVcm" 



Ms 



10i4h^iMo 



2/3 



(5) 



i)'''E(z)-/X-r, 



fg being the baryon fraction. This assumes a metallicity of 
0.3 Zq and a mean molecular weight /i = 0.596. This profile 
appears to fit the results of simulations without radiative 
cooling, with approximately 20% dispersion, especially be- 
yond r > 0.1J?200 (see also MitcheU et al. 2009). 

Assuming the initial gas distribution to have this entropy 
profile, Si{r), and assuming this gas to be in hydrostatic 
equilibrium with a given total mass profile, the initial gas 
pressure profile, Pg,i{r) can be obtained from. 



dPs^ 
dMg 

dr 



GM{^ r) 



47rr2(/i3/5^g^^)r 



Si{r) 



-3/5 



(7) 



where M(^ r) is the total mass inside radius r, given by 
the NFW profile, and Mg{r) is the gas mass inside r. These 
equations can be solved with the boundary conditions that 
Mg = at r = 0, and that the total gas fraction inside the 
virial radius, Mgifvir) /M^i^ = Q.b/^m. 

• (b) Simulations that include radiative cooling, however, 
can also be used to obtain an initial profile that we can 
compare with the observed profile in order to determine the 
additional entropy needed. Loken et al. (2002) found from 
simulations with cooling a universal temperature profile in 
the range 0.04 < r Ir^r < 1, given by. 



r,(r) 



i.33(r) 



(1 + (1.5r/r„,0)'-'^ 



(8) 
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where (T) is the emission weighted temperature of the 
cluster gas. We use an empirically derived scaling relation 
to determine this temperature, from the observations of 
Finoguenov et al. (2001), (although there may be an order 
unity difference between these two definitions of emission 
weighted temperature) 



M500 = (2.641° j^)10" Me 



kb{T) 
IkeV 



(9) 



We can then compute the initial gas density profile assuming 
the gas to be in hydrostatic equilibrium and using the above 
mentioned temperature profile. 



d 



Pg,ikBTi 



GA/« r) 



(10) 



pg,i{r) dr L p-rUp 

The gas fraction within the virial radius is assumed to 
be equal to the cosmic baryon fraction, f7i,/f2m- The 
inferred gas density and temperature profiles can then 
be used to derive the initial entropy profile Si{r) — 

• (c) One can also use a profile with a polytropic equation 
of state, P oc . Li et al. (2010) has recently used this as 
a fiducial profile, with F — 1.2. The normalization comes 
from the assumption of total gas fraction within the virial 
radius to be same as the cosmic gas fraction, as in the above 
models. 

The initial pressure profiles resulting from these three 
possible initial entropy profiles can be compared with the 
observed universal pressure profile of Arnaud et al. (2010). 
Given the observed pressure profile, one can get to density, 
temperature and entropy profile, by assuming hydrostatic 
equilibrium. We determine the gas density profile, Pg,u{r), 
using hydrostatic equilibrium condition and the observed 
universal pressure (eqn[l}. 



dPu{r) 



GM(s; r) 



(11) 



pg,u{r) dr r2 

We note that the observed pressure profile refers only to elec- 
tron pressure and we should multiply the electron pressure 
Pe = n^kTe, by the ratio of pe/p., where p{= 0.59) is the 
mean molecular weight, and pe{~ 1.14) is the mean molecu- 
lar weights per free electron, their values being appropriate 
for a gas with metallicity of 0.3^0.0 

Once the gas density profile is calculated, the corre- 
sponding temperature profile is then readily computed as 
Tu{r) = Pu{r)/{kbpg,u{r)/{pmp)), and then the entropy 
profile, Su{r) = kBTu{r)/{pg,uir)/{fi,mpf^\ 

We show in Figure 1 the three initial entropy profiles 
for a given cluster mass: the case of universal temperature 
profile (dashed line), that of a polytropic equation of state 
with F = 1.2 (dot-dashed line), and the scaled entropy pro- 
file of Voit et al. (2005) (dotted line). The solid curve shows 
the entropy profile corresponding to the observed universal 
pressure profile of Arnaud et al. (2010). As expected, the 
computed initial profile using radiative cooling (dashed line) 
has a higher entropy than the power-law profile expected 
without cooling (dotted line). It should be noted that these 
curves should be used only for Logio[r /r^ir] ^ —1.5, where 



The total pressure, assuming a common temperature. 



. final (observed) 

powei — law 
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. polytropic 
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Figure 1. The entropy profile of ICM gas, for Mvir = lO^'^h'^ 
M0. The solid line shows the entropy profile corresponding to the 
observed universal pressure profile of Arnaud et al. (2010). The 
dotted line shows the scaled power-law entropy profile from Voit 
et al. (2005), the dashed line shows the initial profile computed 
using the universal temperature profile of Lokcn et al. (2002), and 
the dot-dashed line shows the entropy profile corresponding to a 
polytropic gas with F = 1.2. 

the observational data are more secure than in the inner re- 
gions. We can see that the inferred entropy at r ~ 0.1r„ir 
is several times higher than the one expected from purely 
adiabatic case. 

Given these choices for the initial entropy profile, we 
choose the one with power-law entropy, of Voit et al. (2005) , 
for the following reasons. Firstly, the polytropic gas profile 
lacks physical motivation, since it is not clear why the ICM 
gas should behave in a polytropic manner in the absence of 
feedback. Secondly, although the universal temperature pro- 
file from radiative simulations is an appealing choice, it is not 
clear as to what extent it depends on the parameters used 
in the simulations, and also its possible redshift variation is 
unclear. On the contrary, the power-law entropy profile is 
simple and rests on basic assumptions for the behaviour of 
gas in a gravitational potential well. 



3 INJECTED ENTROPY 

Next we compute the additional injected entropy needed to 
reconcile the initial profiles with the observed ones. The ex- 
tra entropy can be expressed in different ways depending on 
the manner in which entropy is injected into the ICM. Pre- 
vious workers have characterized the additional entropy as 
a constant entropy 'fioor' that is independent of the radius. 
For example, Li et al. (2010) have recently used a value 
of Sinj ~ 100-200 keV cm"^ to determine the ICM pres- 
sure profile for their calculation of the SZ effect. This choice 
of characterizing the injected entropy is motivated by the 
curves showed in Figure 1, in which the gas at lower radius 
appears to have an entropy 'floor'. 

There is, however, one problem with this characteriza- 
tion. When one adds entropy to the initial profile, the gas 
shells move outward. In other words, entropy added at one 
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Figure 2. Same as in Figure 1 except that the profiles are plotted 
against the gas fraction Fg instead of the radial coordinate. 



position is carried elsewhere by the movement of the gas, 
and so the comparison of final and initial entropies at a 
given radial position does not remain physically significant. 
As Li et al. (2010) have also recently noted, the injected 
entropy should be characterized by the gas shell in which 
the injection has occurred. In other words, one should com- 
pare Sf{Fg) with Si{Fg), where Fg{r) = Mg[s^ r)/M^ir is 
the gas fraction inside a radius r. Note that this gas frac- 
tion is different from the usual definition of gas fraction in 
the literature, fg = Mg(r) / r); we distinguish these 
two distinct parameters by using two different notations for 
them. 

We show in Figure 2 the corresponding entropy profiles, 
as functions of Fg, for Mvir = 10^*/i~^ M©. It is clear that 
the observed profile (solid line) corresponds to a case where 
the gas fraction within rvir has decreased from the initial 
value Q.b/Q,m ~ 0.16, implying outward movement of gas 
shells. The entropy 'floor' that appears in the entropy profile 
in radial direction exists here only for the very inner shells 
of gas. For most of the gas in the cluster, the (logarithm of) 
final entropy (solid line) is approximately separated from 
the initial profiles by an additive factor, which translates 
to a multiplicative factor for entropy. Viewed in this way, in 
terms of entropy injection in Lagrangian gas shells, we find it 
useful to characterize the entropy injection as a combination 
of an 'entropy floor' and a multiplicative factor for the initial 
profile: S s{Fg) = + r^S^iFg). 

The two terms have different physical implications. 
Note that the entropy [S = kT/n^^^) used here is propor- 
tional to the logarithm of gas entropy. If there had been no 
constant term (Sc), then we could write the excess energy 
per particle as Ae, and the entropy s — Const x In 5, then 
we have Ae = TAs oc oc T{ri— 1). Here, since the tem- 

perature is a weak function of radius, a constant rj would 
imply a near uniform injection of energy. 

On the other hand, if the entropy injection were to be 
only characterized with a constant term, then we would have 
Ae oc TxSc/Si oc n^^^, which would imply that energy feed- 
back works mostly in the inner regions where the gas density 
is highest. Our adopted characterization, with a constant 



term and a multiplicative factor essentially allows provisions 
for both these modes of energy feedback. 

In order to determine the level of entropy injection 
(namely, the value of rj), we solve the following equations 
for hydrostatic equilibrium for the final gas profile: 



dr 

dFg 

dPi 
dFg 



\SfkJ 



-3/5 



A-Kr^firrip \ Sf ki,/ 
GM(g: r) Mvi, 



(12) 



where Fg = Mg{r)/Mv,,. 

The three variables to be determined here are 
ri, Pf{r), Fg{r), and we need three boundary conditions to 
solve these equations. One is an obvious choice: (1) Fg{r = 
0) — 0. (2) Another obvious choice is that P/(-R5oo) = 
PuiRsoo), so that the final pressure profile corresponds to 
the observed universal pressure profile. (3) Lastly we use the 
value of gas fraction at an outer radius, but this requires an 
explanation. 

For the third boundary condition, we use the value of 
^"9(^500) (i-e., Mg(_R5oo)/Mvir) that corresponds to the uni- 
versal pressure profile. We note that the gas fraction at -R500 
as given by the universal pressure profile is different from 
the ratio of ili,/^m ~ 0.167, which implies that a fraction of 
the original ICM gas has been thrown out of the system, if 
one considers the ratio Qt,/flm to be universal initially. We 
show in Figure 3 the gas fraction implied by the universal 
pressure profile, and compare with other observations. Note 
that observations of gas fraction within i?5oo usually refers 
to /g(7?50o) = -F'g(-R50o) X { m"rI„„) )- I" Figure 3 we show 
the values of /g(i?5oo) as given by the universal pressure pro- 
file (solid line) for clusters at z — and compare them with 
the observed values (dotted line), as reported by Sun et al. 
(2009). 

We also show by the dashed line another ratio, that of 
the gas mass within a radius 2Rvir and the total virial mass 



is found to 



of the cluster. This ratio, Fg{2Ryir) — 
be almost a constant for all clusters, and close to the value 
of Q.b/0.m., shown by the dot-dashed line, with an average 
deviation of order ~ 10%. This interesting result suggests 
that a fraction of the ICM gas, originally inside -R„ir and 
with a universal gas fraction of Q.i,/Q.m, gets ejected due 
to the feedback process, but that the gas does not extend 
beyond a radius 2Rvir ■ 

Entropy feedback predominantly affects the gas while 
keeping the dark matter unchanged. In this case, gas mass 
conservation will ensure that Mg{Rvir,i) = Mg{2Rvir, f), 
where the latter term is the gas mass within 2Rvir af- 
ter the feedback processes. Dividing both sides by Mvir, 
and recalling that the initial profile is normalized by 
Mg{Rvir,i) / Mvir ~ nt/Slm, yiclds the condition for the fi- 
nal profile, Mg{2Rvir) / Mvir = ^b/^m- The choice of the 
final radius to be 2Rvir is slightly arbitrary. However, this 
choice is consistent with the findings of Battaglia et al. 
2010, who show that there is significant contribution to 
the SZ flux from beyond the virial radius, and going be- 
yond 2Rvir does not add to the SZ power. This result al- 
lows us to use a boundary condition, that the gas fraction 
Fg{2Rvir) = ^^"a/^^"''' ^ fib/f^m- We also show in Figure 
3 by two thin dashed lines the results for Fg(2.5Rvir) and 
Fg{1.5Rvir), and the total spread is of order 17% of the 
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Figure 3. Gas fraction within iJsoo (/g = Afg(R5oo)/A<f{ii5oo)) 
is plotted as a function of cluster virial mass, for the universal 
pressure profile (solid curves), for z = 0. The dotted line shows 
the observations of Sun et al. (2009) and the dot-dashed line 
shows the value of The topmost (solid) dashed line shows 

the fraction Fg{2Ryir)^ which is the ratio Mg{2R^ir)/Myir, as 
computed from extending the universal pressure profile to 2i?,„i,. . 
Interestingly, this line is close to the Universal gas fraction to 
within 10 — 15% for all masses. The pair of thin dashed lines 
above and below this curve are for Fg(2.5Rvir) and i^g(1.5i?„ir)i 
respectively. 



values of Fg{2Rvir). This is the level of uncertainty in our 
results that stems from the lack of precise knowledge about 
the outer boundary condition. 

We find that with the above mentioned assumptions 
and boundary conditions, the following prescription of en- 
tropy injection works satisfactorily in converting the initial 
entropy profile to the final profile as given by the universal 
pressure profile at 2 = 0, for lO" ^ MvIp/Mq ^ 10^^: 

Mvir 



Si„j{Fg) = 200keVcm"^-f2x( 



10" /i-i Mp 



The values of the coefficients provide the best match 
within an accuracy of ~ 15%. This form of the feedback 
implies that low mass clusters need overall more feedback 
than rich clusters, although (from the constant term) all 
clusters appear to need feedback in their central regions. 

Physically these two modes may refer to preheating and 
AGN/SN feedback processes, one in which all the gas parti- 
cles in the ICM are energized and another, in which particles 
in the central region get more entropy than elsewhere. We 
show below that from the comparison of the final pressure 
profile and the initial pressure profiles, that it would be im- 
possible to achieve the final pressure profile with only a single 
mode of entropy injection. 

We show in Figures 4 and 5 two examples of the re- 
sult of this type of entropy injection, for Al^h- = 10^''h~^ 
Mq and Mvir = 10"/i"^ Mq, respectively, at z=0. The top- 
left panels of both figures show the density profiles: dashed 
lines show the initial profile of power-law entropy; the dot- 
dashed lines show the profile corresponding to the universal 
pressure profile and the dotted lines show the final profile, 
which is a result of the entropy injection according to the 
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Figure 4. Comparison between the initial and final profiles, 
before (dashed) and after (dotted) entropy injection, and those 
obtained from universal pressure profile (dot-dashed), for Mvir = 
10^''h~^ Mq. The top-left panel shows the density profiles, the 
top-right panel shows the temperature profiles, the bottom-left 
panel shows the pressure profiles and the bottom-right panel 
shows the entropy profiles. 



above mentioned prescription. The top- right panels show the 
temperature profile, the bottom-left panels show the pres- 
sure profile and the bottom-right panels show the entropy 
profiles. 

Before discussing the results, let us summarize the 
method of our calculation. We have began with a power-law 
entropy profile, and used a prescription for entropy injec- 
tion, which along with the above mentioned boundary con- 
ditions yield a pressure profile. We iterate a loop until the 
three boundary conditions are matched (within 10%). This 
leads to a match with the observed universal pressure pro- 
file because of the chosen entropy injection prescription (|13p . 
As a result, we find in Figures 4 and 5 that the final pres- 
sure profiles (bottom-left panels, dotted lines) match well 
the observed profile (dot-dashed lines). The match for the 
temperature (top right panels) and entropy (bottom right 
panels) profiles are less satisfactory, but we note that the 
dot-dashed lines here are not observed profiles, but inferred 
from the universal pressure profiles. Since we have decided 
to peg our boundary conditions on the observed pressure 
profiles, we conclude that the above entropy injection pre- 
scription gives a satisfactory match with observations. Since 
a single prescription appears to work for the whole range of 
cluster masses, one could call it a universal entropy injection 
prescription. 



4 EVOLUTION WITH REDSHIFT 

4.1 Evolution of entropy and pressure 

Our exercise has so far been limited to observed clusters 
at low redshifts: we have used a physically motivated ini- 
tial gas profile, and compared it to the observed profile to 
deduce a prescription for entropy injection. It is true that 
the entropy injection prescription allows us to speculate on 
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Figure 5. Same as Figure 4a, but for A/vir = 10"/i" 
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the kind of feedback process that may give rise to the en- 
tropy enhancement. The prescription, however, remains a 
phenomenological one, and we wish to explore its usefulness 
beyond the regime within which it has been deduced in the 
previous section. 

In particular, we wish to study the ICM in high redshift 
galaxy clusters, and determine the amount of feedback that 
may have occurred at high redshift. The properties of the 
ICM at high redshift have been studied through the observa- 
tions of evolution of cluster scaling laws, e.g., the relations 
between X-ray luminosity and temperature, and with the 
help of numerical simulations. In the rest of the paper, we 
show that our semi-analytical formalism can help in this 
regard, by determining the amount of feedback needed to 
produce an ICM profile that is consistent with redshift evo- 
lution of cluster scaling relations. Voit (2005) had suggested 
that cluster scaling relations would differ at high redshifts 
owing to the effects of feedback and radiative cooling (see 
also Muanwong, Kay & Thomas 2006; Short et al. 2010). 

At present, however, the usefulness of this method is 
limited by the lack of accurate determinations of scaling rela- 
tion evolutions. Ettori et al. (2004) found that the evolution 
oi M — T relation in a sample of 28 clusters at 0.4 < z < 1.3 
followed the self-similar model, expected from simple cos- 
mological considerations. But they found the normalization 
of the Lx — T relation to decrease with redshift, which was 
later confirmed by Hilton et al. (2007) from a study of a high 
redshift cluster with XMM — Newton (see also Branchesi 
et al. 2007). Morandi et al. (2007) found negative evolution 
in the M — T relation but a positive evolution for Lx — T 
relation (see also Maughan et al. 2007). At the same time, 
Maughan et al. (2006) found little evidence for any such red- 
shift evolution apart from that expected from self-similarity. 

One reason for this uncertainty is the scatter produced 
by the large X-ray luminosity at the centres of cool core 
clusters (Chen et al. 2007). Cluster mergers also produce 
deviations (Kravtsov et al. 2006). Also, the sample of small 
number of clusters at high redshift suffers from selection 
effects (Pacaud 2007). It is hoped that a large and homo- 
geneous sample will be gleaned from the upcoming XMM 
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Figure 6. The pressure profiles for A/vir = 10^^h~^ Mq at 
z = 0, and 1 are plotted against r/R^oo (in cgs units), for the 
case of constant entropy enhancement. The solid line refers to 
the z = profile, and the thick dashed line is for z = 1. The thin 
dashed fine show the universal pressure profile scaled by a factor 
E{z)^^^ for comparison. Wo also show with dotted line the case 
for the entropy injection scaling as (1 -|- z)~'^'^, which matches 
closely the self-similarly scaled universal pressure profile. 



Cluster Survey (XCS) (Romer et al. 2001), and allow a bet- 
ter determination of the evolution of X-ray scaling relations. 

Keeping these uncertainties in the observed relations in 
mind, we would like to extend our phenomenological study 
to high redshift clusters. For this exercise, we need to as- 
sume a redshift evolution of the feedback prescription. Also, 
instead of changing the boundary conditions with redshift, 
we assume that the boundary conditions at high redshift is 
the same as at 2 = 0, namely, Fg{r = 0) = 0, and that 
Fg{2Rvir) = ^Ib/^m- This will yield the gas profiles, which 
can then be compared with the expectations from scaling 
relation evolutions. This process can then be iterated to de- 
termine the amount of feedback that is required to explain 
the evolution of the cluster scaling relations. This iteration 
will make sense with accurate determination of scaling re- 
lation evolutions, and in the absence of such observations, 
we only demonstrate the first part of this process, namely, 
prediction of ICM at high z corresponding to a given feed- 
back prescription and how this can be used to compare with 
observed cluster scaling relation evolutions. 

We have used simple modifications to the entropy injec- 
tion prescription. Here, one can work out a number of possi- 
bilities, parameterized by a certain set of variables and test 
the resulting pressure profiles against data. In this paper, 
we demonstrate the utility of our approach by calculating 
the cluster scaling relations with redshift resulting from a 
simple redshift scaling of the entropy injection prescription, 
namely, Si„j{z) = (l + z)" Sinj,z=o- This means that we mul- 
tiply the expression of injected entropy in eqn[T3]by a factor 

(i+zr. 

Assuming this feedback prescriptions, we calculate the 
ICM profiles for clusters at high redshift for different values 
of the power-law index a. In Figure 6, we show the result- 
ing pressure profiles for a = 0, for Mvir = 10^^h~^ Mq at 
z — (thick solid line) and 1 (thick dashed line). The thin 
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Figure 7. The evolution of the normalization CI of Lx — T 
relation is shown with redshift, for universal pressure profile 
with self-similar scaling (solid line) and three cases of enhanced 
pressure profile: for a constant entropy injection at all redshifts 
(long dashed line), for Sinj oc (1 + 2:)""'^ (dot-dashed line), and 
Sinj OC + (short dashed line). The three dotted lines show 
the results of simulations GO, FO and PC (see text for details). 

line shows the self-similarly scaled universal pressure pro- 
file at z = 1, scaled by the factor E{z)^^'^ . The resulting 
profiles clearly have larger values for pressure than those 
expected from a self-similar evolution of the pressure pro- 
file. We can similarly calculate the pressure profiles for other 
cluster masses and at different redshifts. 

We also show in Figure 6 the case of a = —0.8, for a 
negative evolution of entropy injection, at z = 1 and we find 
that this matches the self-similarly scaled universal pressure 
profile closely. We have checked this at other redshifts z < 1 
and found that they match within an accuracy of 25%. Our 
calculations with different values of the parameter a show 
that the pressure profiles at high redshift are increasingly 
suppressed with increasing values of a, as expected. This 
will result in different forms of evolution of various cluster 
scaling relations, which can be tested against data, and we 
now turn to a discussion on this issue. 

4.2 Evolution of scaling relations 

How do these profiles compare with the observations of clus- 
ter scaling relations? Following Short et al. (2010), we study 
the normalizations of the — T and the SZ decrement 
Yx — M relation (see below) and their evolution in redshift. 
For the Lx —T relation, we assume a power-law relation, 

Lx E{z)-^ = Cl{Tsi/6keVT , (14) 

where Lx is the X-ray luminosity at 0.5 — 1 keV band and 
Tsi is the spectroscopic like temperature, as defined by Maz- 
zotta et al. (2004). One can study the evolution of the nor- 
malization CI, with redshift, for ICM profiles given by self- 
similarly scaled universal pressure profile, and those with 
other forms of evolution. 

In Figure 7 we show the redshift evolution of CI, in the 
units of 10''^/i~^ erg s~^. The solid line shows the evolu- 



tion obtained from a self-similarly scaled universal pressure 
profile, which means that the profile is scaled by a factor 
E[z)^^^ . The long dashed line shows the case of entropy 
injection being constant with redshift, the dot-dashed line 
shows the case of the scaling (1 -I- z)~^'^ , whereas the short- 
dashed line shows the results of Sinj oc (1 -I- z)"'^. 

This figure can be compared with Figure 12 of Short et 
al. (2010), which shows the results of different simulations 
and observations. As mentioned earlier, the data has large 
scatter and the form of evolution remains unclear. Simula- 
tions with feedback, however, show either positive evolution 
(from FO simulation including AGN and supernova feed- 
back, by Short et al. (2010)) or negative evolution (PC sim- 
ulation, including preheating, by Hartley et al. 2007). Sim- 
ulations without any feedback (GO, in the Millennium Gas 
simulations series; see Short et al. 2010) show a mild positive 
evolution. However, the amplitude without feedback fails to 
match observations. Thus some form of energy injection is 
absolutely necessary. 

We show in Figure 8 the redshift evolution of the nor- 
malization C2 of the relation between Yx — M500 relation, 
normalizing at M = 5 x 10^*/i"^ M©, 

It has been recently shown that Yx = MgT (evaluated within 
a certain over density, here R500) behaves as a low-scatter 
proxy for the cluster mass, regardless of the dynamical state 
of the cluster (e.g., Kravtsov et al. 2006), since it essentially 
expresses the total thermal energy of the gas which is not 
strongly affected by cluster mergers. 

The parameter C2 has the units of keVMQ . It has been 
shown that the universal pressure profile is consistent with 
SZ observations at low redshift (Arnaud et al. 2010). We 
therefore do not show the scaling law at 2 = 0, but focus 
on the evolution of the normalization with redshift. Again, 
the curves in Figure 8 shows the expectations from a self- 
similarly scaled universal pressure profile (solid line) and 
those from entropy injection at high redshift with the follow- 
ing scalings: oc [l + zf'^ (short-dashed line), constant in red- 
shift (long dashed line), and oc {l + z)~"'^ (dot-dashed line). 
The dotted lines again show the results of three types of 
simulations mentioned above, without feedback (GO), with 
AGN/supernova feedback (FO) and with pre-heating (PC), 
as reported in Short et al. (2010). 



5 DISCUSSION 

Figure 7 shows that the self-similar evolution of the uni- 
versal pressure profile (solid line) is consistent with the 
available data. The case for constant entropy injection (dot- 
dashed line) and positive evolution do not fare well, and the 
match with the data improves for negative redshift evolution 
(dashed line) of entropy injection. As mentioned earlier, the 
self-similarly scaled universal pressure profile is equivalent 
to a scaling of entropy injection as (1 -I- z)~°'*, and the solid 
line therefore also shows the case of a steep negative evo- 
lution with a = —0.8. The results from enhanced entropy 
are largely consistent with the range of results from current 
simulations including feedback and preheating (shown by 
the dotted lines marked as FO and PC). 
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Figure 8. Same as Figure 7, but for tlie normalization parameter 
C2 of Ix - A/500 relation. Data points from REXCESS (Pratt et 
al. 2009) are shown as blue and high redshift data from Maughan 
et al. (2006) are shown with red points. 



The curves in figure 8 show that the self-similar evolu- 
tion of the universal pressure profile is consistent with the 
available data, and that the data favour negative evolution 
of entropy injection in general. The enhanced profile with 
our simple scaling law are consistent with the results from 
simulations including feedback. Comparing with the data, 
which include low-redshift observational data from REX- 
CESS (Pratt et al. 2009; Short et al. 2010) and high red- 
shift data from Maughan et al. (2006) , we find that the case 
of constant entropy injection matches better and, as in the 
case of Lx—T relation, the situation improves with negative 
redshift evolution. A more detailed analysis of the likelihood 
of the values of a that is allowed by the data is beyond the 
scope of the present paper, and we will pursue this approach 
in a future paper, and determine the possible evolution of 
entropy injection, parameterized by certain variables whose 
values could be determined by comparing with existing data. 
But the trend shown in Figure 8, indicates that the avail- 
able data require an evolution of the entropy injection that 
scales as (1 -I- z)" with a ^ — 1. Note that owing to the 
non-gravitational nature of the feedback processes, the en- 
hanced entropy need not evolve self-similar ly with redshift. 
However, the entropy of gas as determined by gravity alone 
is expected to evolve self-similarly, namely, 5* oc E{z)~'^^^ for 
clusters with same mass. Nevertheless, we can compare the 
entropy scaling inferred above with the self-similar scaling. 
In A-CDM model, £(2)"^''^ roughly scales as (1 + z)"'^ "* in 
this redshift range. Therefore, the scaling Si„j oc (1 -I- z)^^ 
can be written as Si„j oc S(^)-2/^(l-H^)°■^which is indica- 
tive of a non-standard evolution. 

A more important conclusion that we can draw from 
the previous section, however, is that that in the presence 
of accurately determined cluster scaling relation evolution, 
as is expected in the near future, one will be able to com- 
pare the predictions from entropy injection with that from 
redshift evolution of cluster scaling relations, and therefore 
iterate in order to determine the entropy enhancement at 
high redshift. 



Fang & Haiman (2008) have argued, on the basis of a 
simple model of evolving entropy floors for the ICM with 
redshift and comparing the results with observed evolution 
of scaling relations, that the entropy floor at high redshift 
is likely to be shallower than that at present. They argued 
that a deceasing entropy floor with redshift is expected from 
feedback processes due to active galactic nuclei (AGN) and 
supernovae driven winds. 

It is however difficult to interpret the amount of feed- 
back only from the level of the entropy floor, since the de- 
fault entropy at high redshift is also likely to be small, from 
cosmological evolution. Entropy, as deflned by S = T/n^^^ , 
would normally evolve as E{z)~^^'^, for clusters with same 
virial mass, in the absence of any feedback process. Any 
deviation from this can be identifled as the result of a feed- 
back process. For example, Ettori et al. (2004) found that for 
clusters with same temperatures, the entropy floor evolves 
as E{z)*^^ S' oc (1 -I- z)°-^. Nath (2004) interpreted this as 
the result of increasing feedback energy deposition at higher 
redshifts, possibly as a result of increase in the abundance 
of AGNs at high redshift. 

Our work in this paper presents a method of determin- 
ing the amount of feedback and its evolution with redshift, 
by using an initial entropy profile and injecting entropy into 
the profile (defined by radially distributed gas shells). It re- 
mains to be seen how the entropy injection prescription is 
related to any physical process of feedback, and we will re- 
turn to this issue in a future paper. 

Finally, we note that our results are accurate up to 
~ 15%, given the accuracy with which we have implemented 
the boundary conditions in arriving at our results. Our anal- 
ysis presented here has ignored the contribution from non- 
thermal pressure of the gas, which can lead to underesti- 
mating the mass by 15 — 20% around r25oo at ~ 25% at the 
virial radius (see for example, Rasia et al. 2004), but which 
would have introduced additional free parameters, and more 
uncertainties in the interpretation of cluster scaling relation 
data, the results of the entropy injection prescription. The 
simplicity of the entropy injection prescription presented in 
this paper, however, may make it possible to include non- 
thermal pressure. 

There is another source of uncertainty that arises from 
the fact that a fraction of the baryons in clusters would 
form stars, and would thus affect the boundary condition 
for normalization. Although Gonzalez et al. (2007) found 
that the stellar mass fraction in a cluster of M500 ~ 10^* 
M0 can be as large as 0.25Qb/^m (see also Giodini et al. 
2009; Andreon et al. 2010), Balogh et al. (2008) have found 
that such a large fraction, if confirmed, would cause prob- 
lem for the paradigm of hierarchical structure formation, let 
alone being an 'overcooling' problem. On the other hand a 
conservative estimate from Bower et al. (2006) is that the 
stellar fraction is ~ 0.1 times the cosmic baryon fraction. 
This is also in agreement with hydro simulations of clusters 
(Ettori et al. 2006). Recently, Simionescu et al. (2011) have 
found from Suzaku observations of the Perseus cluster that 
the baryon fraction approaches the Universal value much 
earlier than the virial radius. All of these would imply that 
our boundary conditions used here would be modified by a 
factor ~ (1 — 0.1), and the resulting profiles would differ 
by ~ 10% and, hence, would decrease the required entropy 
injection by a similar fraction. We have neglected the contri- 
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bution of star formation in this paper because our entropy 
injection prescriptions are accurate up to ~ 15%. 

We would like to thank the referee for a very careful 
reading of the manuscript that has helped us in improving 
the paper. We thank Satej Khedekar for help with Figures 
7 and 8. 
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